Super-Hubble Nonlinear Perturbations During Inflation 
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We show that the non-linear evolution of long-wavelength perturbations may be important in a wide 
class of inflationary scenarios. We develop a solution for the evolution of such nonlinear perturbations 
which is exact to first order in a gradient expansion. As a first application, we demonstrate that 
in single field models of inflation there can be no parametric amplification of super-Hubble modes 
during reheating. We consider the implications of the solution for recent discussions of the back- 
reaction effect of long wavelength perturbations on the background geometry, give a new derivation 
of the equation of motion of stochastic inflation, and demonstrate that if the (generalized) slow- 
rolling condition is not satisfied, then inevitably long wavelength vector modes for gravitational 
fluctuations will be generated. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The hypothesis that causally unrelated regions of 
space-time evolve independently is one of the corner- 
stones in the development of relativity and one of the 
basic assumptions of all relativistic theories. As an impli- 
cation of this hypothesis in cosmology we expect that per- 
turbations with long wavelengths^] evolve independently 
at different points of space, as if they were independent 
homogeneous patches of the Universe (see e.g. |IJ and ref- 
erences therein). This expectation has been explored in 
the linear theory of cosmological perturbations. In this 
context, it is possible to derive integral constraints (see 
e.g. ||) which imply that in the case of adiabatic fluc- 
tuations there can be no generation of fluctuation modes 
on super-Hubble scales. However, it is possible for en- 
tropy fluctuations to be generated on such scales, one 
well known example being the formation of cosmic defects 
during a phase transition which leads to super-Hubble 
structures (see e.g. Another example is the gener- 
ation of axion fluctuations (see e.g. Q). Recently the 
super-Hubble range effects of entropy fluctuations have 
been widely studied in the context of inflationary reheat- 
ing (see e.g. [|-§). 

The aim of this work is to investigate the process of 
structure formation in the non-linear regime, for infla- 
tionary models based on a single scalar field. In Section 
2, we show that the presence of nonlinearities is inevitable 
in a wide class of inflationary models. In Section 3 we find 
a solution (which is exact to leading order in the gradi- 
ent expansion) for nonlinear modes in the case when (but 



not only when) perturbations satisfy a generalized slow- 
rolling condition. This condition is naturally satisfied 
during inflation. For our solution we can define a gen- 
eralized Bardeen parameter which is constant on scales 
larger than the Hubble radius. As a first application 
(Section 4), we demonstrate that there is no amplifica- 
tion of long wavelength gravitational fluctuations during 
reheating. In Section 5 we discuss some implications of 
our solution for the back-reaction of perturbations on 
the background geometry. In Section 6, the Langevin 
equation for stochastic inflation Js) is re-derived in a way 
which includes the effects of gravitational (and not just 
matter) fluctuations, using the results of Section 3 which 
can be self-consistently applied in an inflating Universe 
with nonlinear super-Hubble modes. Finally in Section 
7 we demonstrate that in the case in which the general- 
ized slow-rolling condition is violated, the generation of 
vector modes is inevitable. 



II. GENERATION OF FLUCTUATIONS 

Perturbations are generated in the de-Sitter phase of 
an inflationary Universe as a consequence of quantum 
vacuum fluctuations (see e.g. || for a comprehensive re- 
view). The analysis is based on the consistent quanti- 
zation of the linearized metric and matter perturbations 
in a classical expanding background space-time. In this 
framework, it can be shown that the Sasaki-Mukhanov 
parameter v defined as 'h 



v = a [S<p + = ^C, 



(1) 



*Throughout this paper, we use the terms super-Hubble, 
infrared and long wavelength interchangeably. All refer to 
modes with wavelengths larger than the Hubble radius. 



^The cosmological scale factor is denoted by a(t), the Hubble 
expansion rate by H , and the scalar matter field by ip. 
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in the linear regime obeys the equation of motion of a 
free scalar field. Here, C is the Bardeen parameter Jl2| |: 



(2) 



where <fi is the variable which describes scalar metric per- 
turbations in longitudinal gauge in a spatially flat Uni- 
verse: 



where the subscript IR indicates that we only include 
infrared modes. 

Integrating the equation of motion for the background 
field and the Friedmann equations in the slow rolling 
regime gives: 



*Vdtp 3 k 
—— = m p log{ — 

v P v 



ani,, 



(7) 



ds 2 = (1 + 2(j))di 2 - a 2 (t)(l - 2(f))5 l jdx l dx i , (3) where m p = (8irG) i is the Planck mass and V(ip) is the 



The scalar sector is the dominant one for cosmological 
perturbations generated through slow-roll inflation (see 
e.g. |TJj and references therein). 

The action for linearized metric and matter perturba- 
tions about a classical expanding background becomes 
the action of a free scalar field v with time-dependent 
mass (the time dependence is determined by the back- 
ground geometry and is negligible on scales smaller than 
the Hubble radius). Hence, the perturbations can be 
canonically quantized, and expectation values of v (and 
hence also of £) can be easily calculated (see e.g. ||) once 
the state of the system is specified. The state can be spec- 
ified mode by mode in Fourier space. Our choice is the 
following: We fix the initial conditions for each Fourier 
mode when its wavelength is equal to the Planck length 
[[], and assume that we have vacuum expectation values 
at that time. This assumption is very similar to the usual 
choice of a vacuum state at the onset of inflation, since 
- at least for the usual dispersion relations @,^5| - the 
wavefunction is almost constant on sub-Hubble scales. 
The advantage of our choice compared to the usual one is 
that we do not extrapolate the physics above the Planck 
energy. However, for the usual dispersion relations for a 
free field, the difference is immaterial. Using the equa- 
tion of motion for £ in the inflationary era, we obtain the 
following expression for < £ 2 >: 



< C > 



d 3 k 



Hp 



(2tt) 3 2fc 3 ¥>p 



a 2 H 2 



(4) 



Note that C contains the full information about the lin- 
earized metric (and matter) fluctuations (valid as long as 
the slow-roll approximation holds). 

Now we turn our attention to the infrared part of the 
spectrum <C 1) for which 



-iL c . 



Thus 



< 



>IR- 



IR 



d 3 k H A p 
J2^f2k^ 



(5) 



(6) 



*The values of the variables at this time are indicated by a 
subscript p, and they are functions of k. 



scalar field potential. Assuming a power law potential: 



we get: 



<ff I = tp 2 - 2m 2 alog( ) 

y " am„ 



(8) 



(9) 



Substituting in (6) and carrying out the integral gives: 



< 



>IR- 



a 



48(4 + a)7r 2 m p m 



(^)4 ( ^ n ^)4+ Q; (10) 



for | ip <C \<fii\, where ipi is the background field value 
when inflation starts. 

We see that if inflation lasts long enough, i.e. if fa is 
large enough, then < 4> 2 >ir can become larger than one 
and thus infrared perturbations can go nonlinear. As a 
matter of fact, although the amplitude of each mode may 
be small, due to the large phase space of infrared modes, 
the actual value of 4> i n real space may be greater than 
1, and so we may not be allowed to use perturbation 
theory to expand the Einstein equations. The possible 
importance of nonlinear effects for the evolution of space- 
time was pointed out in [p^UT^ I, though with a different 
interpretation. We will return to this issue in Section 5. 

Let us now study under which conditions < 4> 2 >ir is 
indeed nonlinear. As a first step, note that the value of M 
is constrained by the observational data. The amplitude 
A of a linear cosmological perturbation mode (measured 
in terms of £) when it enters the Hubble radius after 
inflation is the same as when it exits the Hubble radius 
during inflation ^|. Consider, for example, the largest 
wavelength observable mode, the mode which is entering 
the Hubble radius today and is given by: 



k = a H a ~ ( 



M 



Tcmb m 



)[—)[a e m p ), 



(11) 



§ Although this statement is conventionally taken to be true, 
there are exceptions as has recently been shown for a cer- 
tain class of inflationary models (see Section 4 and references 
quoted there). 
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where the subscript refers to the value of quantities 
today, Tcmb denotes the present temperature of the mi- 
crowave background, and a e is the value of the scale fac- 
tor when inflation ends. Then, combining (4) and (9) 
gives us the amplitude of the mode: 



A = 



Hp 

Vp 2 



M 



(— nM 



TcMBfn, 



MH Q 



Thus we end up with: 



M .1, 1, „, °+2 

-^AH67--lo 9 A]-^. 



(12) 



(13) 



The value of A is constrained by various observations, 
in particular by the amplitude of the cosmic microwave 
anisotropies on large angular scales Il8| , which give A ~ 

io- 10 . 

A theoretical requirement for a successful inflationary 
scenario is that all the observable modes must have exited 
the Planck length during inflation. The consequence is 
that (pi is larger than <p p for all the observable modes. 
Using (9), (12) and (13), this requirement becomes: 



ipi > m p [log( 



MH ' 



m p [67- -log A}?, (14) 



(for this order of magnitude estimate we can set a = 1). 
Inserting this constraint in (10) gives: 



< <f > IR > A[67 - \logA] z+ % 
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-4. 41+0. 93a 



. (15) 



at the end of inflation, for A ~ 10~ 10 . We see that 
even if inflation starts at the latest possible time which 
is allowed by theoretical/observational considerations (it 
is usually taken to start much earlier), < 4> 2 >ir almost 
reaches the nonlinear regime for a massless theory (i.e. 
for a — A) just before reheating starts. If inflation starts 
at values of cj) for which the energy density is comparable 
to the Planck density, then the phase space of infrared 
modes is much larger and the nonlinear regime will be 
reached much earlier. 



III. GRADIENT EXPANSION 



We use the following ansatz for our metric: 



ds 2 



^dt 2 



(16) 



which is a generalization of a metric with linear scalar 
perturbations. This ansatz, though not a general solu- 
tion, represents the most important sector of the metric 
since any such metric can be continuously connected to 
a metric with linear scalar perturbations which is the 
dominant sector in a Universe with linear perturbations 
generated during slow-roll inflation jl3). We also con- 
sider the simplest model of matter which has one scalar 
field (the inflaton). 



The Lagrangian and energy-momentum tensor are 
given by: 



C=^d a cpd a cp-V(<p), 
T£ = d»(pd u <p-£6Z. 



(17) 
(18) 



In writing the Einstein equations we keep the leading 
terms in their gradient expansions, assuming that the 
wavelengths are large. Naturally, in the Goo and Gy 
equations, all the gradient terms drop out Q and we are 
left with the Friedmann equations for the above metric: 



2 8irG.(p" 



-(- 



V(<p)), 



-2ip" 



3 v 2 

3^ 2 = 8^G(-^- + U(^)), 



(19) 
(20) 



and of course the more interesting Goi equation: 

{ip'),i = AxG<pf(p ti (21) 

where ' denotes e - ^-^-. 

We immediately sec that not every ip can satisfy the 
last equation since the right hand side must be a perfect 
gradient. The condition that tp must satisfy is: 



(<p'),i<P,i = 
or, equivalentlyp|: 



(22) 



(23) 



From Equation (23), we can see that surfaces of con- 
stant ip' and constant cp must be tangential to each other 
and consequently are the same. This means that tp' must 
be a function of ip, at least in the finite regions of space 
in which the gradient of (p does not vanish: 



(24) 



where g(tp, t) can be any function of ip and t. Then the 
Goi equation reduces to PI 



<// = ^Gg{p,t). 



(25) 



**It is possible to show that the Gij equations for i j can 
be satisfied only if we include the vector and tensor sectors 
of the metric perturbations. However their magnitudes go to 
zero in the long wavelength limit and hence, in this limit, they 
do not affect the equations for scalar perturbations. 

^Equation (22) implies that the cross product of two gra- 
dient vectors vanishes, so the vectors must be in the same 
direction. 

The integration constant can be absorbed in g((p,t). 
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Since there is no explicit time dependence in this formal- 
ism, we assume that g has also no explicit time depen- 
dence: 



l(<P,i) = g{tp). 



(26) 



From now on, we call Equation (24) (with no explicit 
time dependence) the generalized slow-rolling condition 
due to its resemblance to the field equation during the 
slow-rolling phase. Note that at this point, we have not 
shown that the absence of time dependence is required 
in (p6|), but simply that this ansatz gives a consistent 
solution. However, below we show that in the case of 
single field inflationary models which initially satisfy the 
slow rolling condition, (|26| ) is indeed satisfied. 

As one can easily see, in the linear regime Equations 
(25) and (26) lead to the vanishing of the non-adiabatic 
pressure defined as: 



Spnad = ^ 



p p 



(27) 



This is the condition derived in 
£ for the perturbations. 

Combining Equations (24) and (26) gives: 

ip" = AwGip' 2 , 



for the constancy of 



(28) 



which can also be obtained by combining the Friedmann 
equations (19) and (20). 

After neglecting the gradient terms and using Equa- 



tions (24) and (26), both the G 00 
duce to: 

4vrG g' 2 - 3(4ttG g) 2 



and Gij equations re- 



~8ttG V(tp). 



(29) 



Thus we see that, in the long wavelength limit, the 
Goi equation is consistent with the Friedmann equations 
given that the field tp satisfies the generalized slow-rolling 
condition (24). 

As a matter of fact, this is what happens during slow- 
roll inflation, when we have: 



1 dVjip) _ 1 dV(ip) 
■iH dip ~ Zip' dip 



(30) 



and hence (24) is satisfied. 

Because of this property, even after the end of infla- 
tion, as long as the perturbation is in the super-Hubble 
regime, tp' remains a function of tp. The reason is that, 
using the field equation, we can track the evolution of tp 
and tp' into the reheating era, as functions of the initial 
ip during inflation. By eliminating the initial tp, it is pos- 
sible to find ip' as a function of tp at any time. However 
it will not be a single-valued function when tp is oscillat- 
ing. Figure ([!]) shows the behavior of g(tp) (solid curve) 
during the reheating era, obtained by solving (29). The 
dashed curve is y/V((p)/(6nG), the asymptotic limit of 
g{tp) for tp » m p . 




FIG. 1. The solid curve is a solution for g(tp) during the 
reheating era and the dashed curve is y/V(tp)/(6irG). The 
horizontal axis is in units of Planck mass while the vertical 
axis has arbitrary units. 

Dividing Equation (24) by Equation (25) and integrat- 
ing the result with respect to tp, we find that: 



C(x) = ip- 4ttG 



(31) 



is an integration constant and as such does not depend 
on time. 

It can easily be verified that for linear perturbations 
8£(x) reduces to the usual Bardeen parameter £, which 
is known to remain constant for adiabatic super-Hubble 
modes ^} 



*C(aO = Sip- 4ttG^ 
9 



Sip 



g5{e-H) 



Sip - ^(H5cf> + Sip) = C(x). 
H 



(32) 



Thus, we have been able to generalize the Bardeen pa- 
rameter to the case of nonlinear long-wavelength fluctu- 
ations. 

Let us now consider the case of multiple scalar fields. 
In this case, Equation (21) changes to: 



{ip') ti = 4acG<ff a <p a . 



(33) 



where the subscript a denotes the field index and is 
summed over, if repeated. The consistency condition (22) 
takes the form: 



^In linear perturbation theory, S<j> is forced to be equal to Sip 
by the off-diagonal Gij equations for all theories in which the 
off-diagonal components of Tjj vanish to linear order, which 
is in particular the case for scalar field matter. 
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{f' a ),i<Pa,j = {<f'a) ,jVa,i- (34) 

We see that if the fields satisfy the generalized slow rolling 
condition 

V' a = J~9(<P,t), (35) 

then (34) is satisfied and we again arrive at (25). How- 
ever, for the multiple field problem this result is not gen- 
eral. Note that (35) is only a sufficient condition. 

During inflation, if the effective mass of a field compo- 
nent is larger than the Hubble constant, it is damped in 
a few e-foldings. Otherwise, it slow-rolls down the poten- 
tial. So a few e-foldings after the onset of inflation, all 
the remaining field components undergo slow-rolling and 
satisfy (35). Wc conclude that during multiple- field in- 
flation, we can still use the above results. However when 
the inflation ends, (35) breaks down and, unlike the sin- 
gle field case, we end up with a mixture of scalar and 
vector metric perturbations, as described in Section 7. 

We are also able to define the generalized Bardeen pa- 
rameter (31) for each slow-rolling path in field space, 
since the evolution is effectively of single field nature 
along each path. However, the parameter does not re- 
main constant after inflation ends since (35) breaks down. 

In conclusion, we have shown that causally discon- 
nected regions of the Universe evolve independently in 
the sense that locally observable parameters satisfy the 
conservation law (31), which is a generalization of the 
conservation of the Bardeen parameter in linear pertur- 
bation theory. As mentioned, our result is only true in 
the case of (in general) non-linear scalar perturbations 
satisfying the slow-roll condition. 

In the next two sections we consider some of the ap- 
plications of this result. 

IV. PARAMETRIC RESONANCE OF 
SUPER-HUBBLE MODES DURING REHEATING 

It was recently suggested [[l9) that parametric reso- 
nance during the reheating phase of an inflationary Uni- 
verse may lead to an exponential amplification of 
super-Hubble scale gravitational perturbations. If true, 
this would affect the usual predictions of inflationary 
model for observables such as the matter power spectrum 
and the spectrum of cosmic microwave anisotropics. 

In Ref. jSlJ it was shown that, although there is no 
causality constraint which prohibits the amplification of 
super-Hubble ( but sub-horizon) modes during reheating, 
the effect does not occur in a simple massive scalar field 
model of chaotic inflation (i.e. with a = 2 ). This result 
was shown to be true even beyond the linear analysis, 
using numerical methods fl22fl . Recently, a general no-go 
theorem for resonance of long wavelength scalar gravita- 
tional fluctuations in the context of a single scalar matter 
field theory was suggested in |23| . In this last reference, 



the effect was also investigated numerically for a matter 
theory with both quadratic and quartic terms (see also 

The problem with the analysis of j23j or with every 
other analytical approach which is based on consider- 
ing the evolution of £ (in the linear regime) through the 
turning point of tp is that C is ill-defined when tp van- 
ishes. However, £ is continuous at this point since the 
integrand of (30) diverges as (tp — tpo)~? ,|***| so the inte- 
gral is continuous. Therefore, £ is well defined through 
the turning. We now argue that £ is in fact constant 
throughout. For the argument we refer back to Figure 1 
which shows a sketch of the time evolution of g(tp) during 
the period of oscillation of tp. During each time interval 
between the turning points, £ and hence £ are constant. 
At the turning points, £ is continuous, and hence it is not 
possible that at these points C undergoes a jump. The 
fact that £ is constant implies the absence of parametric 
amplification of super-Hubble modes. 

In contrast to the single field pointed out at the 

end of Section 3, this formalism cannot be applied to the 
multiple-field case during reheating, and thus we may ex- 
pect non-trivial behavior of the fields during this era. In 
fact, such behavior has been observed for a specific two- 
field potential P,|6| p5|j2r| , where parametric amplification 
of isocurvature modes can lead to non-linear amplitudes 
of infrared modes. We will return to this effect in Section 
7. 



V. BACK-REACTION OF ADIABATIC 
INFRARED MODES 

It has been claimed that the growth of < 4> 2 > 

during slow-roll inflation can lead to significant correc- 
tions to the background Friedmann equations at second 
order in perturbation theory (see also [ p7|j2^1 for sim- 
ilar discussions based on the back-reaction of infrared 
gravitational waves). We claim that, as a result of the 
analysis of Section 3, this effect, though formally valid, 
cannot be identified by local observers as an effect of in- 
homogeneities. The reason is that, as we showed, for 
perturbations generated during slow-roll inflation, the 
local Friedmann equations are always satisfied. Note 
that our analysis includes the effects of the leading back- 
reaction terms in the abovementioned references. Obvi- 
ously, back-reaction terms which are of higher order in 
gradients are not included in our analysis. 

Due to the fact that inflation is followed by a radi- 
ation dominated era, modes which enter the Hubble ra- 
dius after inflation are damped by the Hubble expansion. 



***This can be obtained by considering Equation (21) close 
to the turning point. 
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Consequently, the phase space of constant amplitude per- 
turbations is shrinking after inflation. This effect can 
smooth out the density field so that it becomes of linear 
order again at some point in the radiation dominated era. 
This will happen in models in which inflation does not 
last much longer than the minimal number of e-foldings 
required to solve the problems of standard big bang cos- 
mology. However, in many models of chaotic inflation, 
inflation lasts long enough, so that the remaining phase 
space at the time of equal matter and radiation is large 
enough to give nonlinearity. However, even in this case 
there has not been enough time for the non-linearities to 
enter the observable region (sub-Hubble scales). 

This does not mean, however, that the presence of non- 
linearities has no effect at all. As we will see in the next 
section, the process of generation of perturbations dur- 
ing inflation is influenced by nonlinear effects and thus 
the spectrum is distorted compared to what would be 
predicted in linear theory 



< Cc (xi,*i)Cc (x 2 ,i 2 ) > 
5H 6 



(39) 



w F(Ha\x 1 -x 2 \,H\t 1 -t 2 \ 



where F(0, 0) = 1 and F((3. r) — > as (3 or r goes to 
infinity (more precisely, F(J3,t) <C 1 if /3 3> 1 or r 3> 1)- 
The explicit form of F is complicated and we do not write 
it down here. 

Taking the time derivative of the definition of C (see 
equation (31)), and using the fact that ijj' = —H ~ 



^Y?-V((p) during inflation, we get 

p „ A w d lo g(ffK -i , 



(40) 



We can neglect g' in (29) to find g as a function of V(<p) 
in the slow-roll approximation. Plugging this into the 
above formula, we obtain the coarse-grained equation of 
motion: 



VI. EQUATION OF MOTION FOR STOCHASTIC 
INFLATION 

In this section we will use coarse graining on the scale 
of the Hubble radius to provide a new derivation of the 
equation of motion for stochastic inflation Q which takes 
into account the gravitational fluctuations. 

Let us define the coarse-grained Bardeen parameter, 
£ c; in Fourier space as follows: 



C c (fc) = C(k)W(^) 



(36) 



where the window function W(k) can be chosen to be^ 



h h 2 
aH a z H z 



(37) 



which is smooth on the Hubble scale. In a patch larger 
than the Hubble radius but small enough so the non- 
linearity in £ is negligible, we can use £ (usual Bardeen 
parameter) instead of <5£. Making use of the fact that 
background £ is constant in this patch we have: 



( = sc = (. 



(38) 



We can use our knowledge of the linear quantum gen- 
eration of perturbations (Section 2) to find the statistical 
properties of with the result 



"'This is only one example of a cutoff function. In general, 
one can use any function which is close to one for k < aH and 
tends to zero as k goes to infinity fast enough to eliminate the 
ultraviolet divergence. 



2 u d\og(V) 



. 5 H*, . 

£(x,i) . 

48 tt ' ; 



(41) 



Here, all quantities are coarse-grained, and the variable 
£ is a random Gaussian field whose two-point correlation 
function is given by: 

< £(xi,iiK(x 2 ,< 2 ) >= Ffffolxx-xal.iTltx-ial). 

(42) 

Note that since £ is a metric perturbation variable, 
the above equations correspond to metric back-reaction 
yielding stochastic dynamics. 

Let us compare this calculation with the original one 
in ||. There, instead of the window function (37), a 
step function was used. More importantly, gravitational 
perturbations were neglected. As a result, the numerical 
factors are different. In addition, in the original analysis 
the random variable had a white noise spectrum whereas 
our noise has a finite correlation time. In spite of these 
differences, we expect that the qualitative behavior of the 
solutions will stay the same. 

One major difference between our analysis and the 
usual analysis of stochastic inflation (see e.g. (^9|) is that 
if inflation stops at some point, and consequently the 
generalized slow rolling condition breaks down (see the 
end of Section 3), then, as elaborated in Section 7, the 
evolution takes a nonlocal form. In particular, the end of 
inflation in some region of space can affect the inflating 
regions in an acausal way. 

Limiting the analysis to the linear regime leads to the 
standard results of the stochastic inflationary scenario. 
Namely, since in the large picture different points in space 
undergo independent random evolution, we end up with 
a scale invariant spectrum for £ (or £) with a logarithmic 
correction in Fourier space due to the finite correlation 
length of the random field £ . 
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However, as argued in Section 2, there is a large class 
of scenarios which allow infrared nonlinear perturbations. 
As a matter of fact, the main result of Section 2 is that 
the cumulative effect of the second term in (41) may be- 
come important even when its magnitude is negligible. 
Hence, even when the self-reproduction of the stochastic 
scenario does not take place, infrared nonlincarities may 
be important 

This nonlinearity affects the generation of perturba- 
tions by changing the background parameters which ap- 
pear in the amplitudes as well as the Gaussianity of the 
inhomogeneities. We postpone the calculation of this cor- 
rection to a future work. 



VII. VECTOR MODES 

In this section we investigate the case in which the gen- 
eralized slow-rolling condition (24) or (35) is not satisfied. 
For the general metric: 



ds 2 = N 2 dt 2 -'jijdx'dxi, 
the Oi Einstein tensor elements are given by 



where 



G° = A ':A., A/ :/ ). 



K = K\. 
For the metric used in Section 3: 



N 



mi 



and thus: 



Ki 



(43) 
(44) 

(45) 
(46) 

(47) 

(48) 



The G\ equation is not satisfied in general. It is well- 
known that beyond linear order in perturbation theory, 
tensor and scalar perturbations mix. However, in our 
case we cannot satisfy the equation by adding tensor per- 
turbations. We thus include vector metric perturbations 
in order to satisfy the G\ equation. 

The presence of vector modes would modify K as fol- 
lows: 



(49) 



where Ai is some vector field made up from the metric 
variables representing vector perturbations. Substituting 
this into (44) leads to the following G° equation: 



2Vi// - 2R.A + V x (V x A) = BnGtp'Vip (50) 

where R is the Ricci tensor. If R vanishes, it is easy to 
find A and tp 1 in Fourier space: 



87rG 

A k = — x (k x C k ), 

, -4«ttG 
= — p— (k-C k ), 



where 



(51) 
(52) 

(53) 




C k = 8ttG J d 3 xe" 4k x v3'V^, 

is the Fourier transform of the right hand side of 

One can take the R term to the right hand side of (|5C 
and write a perturbative expansion in R for A and tp' . 
We conclude that for general nonlinear scalar field pertur- 
bations the metric which satisfies the Einstein equations, 
in general, must include nonlinear vector perturbations. 
One known example of this phenomenon is the formation 
of vector perturbations in the case of topological defect 
formation during a phase transition, where super-Hubblc 
structures form in a short time || . Another example con- 
sidered in flsppstl (see the end of Section 4) occurs in in- 
flationary reheating, since during multiple field reheating 
the generalized slow-rolling condition (35) breaks down. 
Note that in the back-reaction calculation of |2^] and in 
models in which the infrared modes lead to nonlineari- 
ties, it is not enough to include only scalar perturbations 
since vector perturbations are generated as well. 



VIII. CONCLUSION 

We have demonstrated that in many models of infla- 
tion, the large phase space of infrared modes leads to 
nonlinearities, even when the amplitude of each Fourier 
mode is small. We were able to find a solution of the Ein- 
stein field equations for nonlinear fluctuations which is 
exact to leading order in the gradient expansion, and thus 
will be accurate to describe the infrared modes. For this 
solution, we were able to define a generalized Bardeen 
parameter which is conserved in time. 

As a first application of this formalism, we were able 
to show that in models with a single matter field, there 
can be no parametric amplification of super-Hubble cos- 
mological fluctuations during inflationary reheating. Ap- 
plied to the problem of back-reaction of infrared modes, 
our solution implies that the back-reaction effect is not lo- 
cally identifiable as an effect due to inhomogeneities since 
the local Friedmann equations are satisfied. We were able 
to use our formalism to give a re-derivation of the equa- 
tion of motion for stochastic inflation which takes into 
account the effects of gravitational fluctuations. Finally, 
we have shown that nonlinearities inevitably lead to the 
generation of vector perturbations if a generalized slow- 
rolling condition is not satisfied. 

This work has been supported in part by the 
US Department of Energy under Contract DE-FG02- 
91ER40688, Task A. 
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